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Abstract 

Ferromagnetism is predicted in undoped diluted magnetic semiconductors illuminated by in- 
tense sub-bandgap laser radiation . The mechanism for photo- induced ferromagnetism is coherence 
between conduction and valence bands induced by the light which leads to an optical exchange 
interaction. The ferromagnetic critical temperature Tq depends both on the properties of the 
material and on the frequency and intensity of the laser and could be above IK. 

PACS numbers: 
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Information processing in electronic devices is based on control of charge flow in semicon- 
ductor materials, whereas non-volatile information storage exploits ferromagnetism, spon- 
taneous alignment of the spins of many electrons. Spintronics^ aims to achieve a merger of 
these technologies, motivating interest in new ferromagnetic semiconductors like (Ga,Mn)As 
and other (III,Mn)V compounds^, and in diluted magnetic semiconductors (DMS) like 
(Cd,Mn)Te and other (II,Mn)VI materials^ with free moments that can be aligned by exter- 
nal fields. The optical properties of these semiconductors are particularly interesting in this 
respect, because laser radiation can control both charge and spin dynamics^. The possibility 
of optically controlling semiconductors with dilute magnetic elements is now being explored 
actively^. 

In this paper we predict a qualitatively new effect, photo-induced ferromagnetism, in 
which magnetic order is induced in an otherwise paramagnetic (II,Mn)VI DMS by laser light 
that is below the absorption edge. The influence of the laser field on the electronic system 
is reactive, rather than dissipative and follows ultimately from a change in the effective 
electronic Hamiltonian analogous to those that yield dipole forces in atomic physics. The 
II- VI parent compounds are intrinsic semiconductors with states of p character at the top 
of the valence band and of s character at the bottom of the conduction band^. They have a 
sizable interband optical matrix element, dcv = (f |er|c). In (II,Mn)VI alloys, Mn substitutes 
for the group-II atoms. The electronic structure of the external s,p shells of Mn and group 
II atoms is very similar, so that the conduction and valence band are barely affected by 
moderate Mn doping. However, the d shell of each Mn atom is only half full, in contrast to 
the full shells of the atoms for which they substitute, and forms a local moment with spin 

= |. These magnetic moments are the only low energy electronic degrees of freedom in 
ideal (II,Mn)VI materials. 

The magnetic moment of a Mn d shell interacts with the spins of electrons in both 
conduction and valence bands through the exchange coupling 

'^exch = ^ JbMi ■ Sb{Ri), (1) 

i,b 

where Mj refers to the Mn atom located at Ri, the spin density of the b = e conduction elec- 
trons and b = h valence band holes is 5b (-Rj), and Jb are the exchange coupling constants^. 
When the Mn atoms are spin-polarized this interaction leads, in a mean field and virtual 
crystal approximation (MFVCA) whose approximate validity is well established^!^, to an 
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effective magnetic field gfisB^^ = Je^hCun-M. experienced by the spins of the carriers, where 
CMn is the density of Mn atoms and Ai is their average magnetization. This leads to band 
spin-splittings as large as 100 meV, evident in optical spectra and known as giant Zeeman 
splitting^. 

In doped systems band electrons mediate RKKY interactions between Mn spins. On 
the other hand, virtual fluctuations in the Mn valence lead to short range superexchange 
antiferromagnetic interactions which couple only nearest neighbors. In a sample with a 
fraction x of Mn randomly located at the cation sites, Xes = x{l — x)^^ is the fraction of 
sites containing Mn atoms without a magnetic first neighbor. Consequently, samples with 
a low Mn concentration (x ~ 0.01) and no band carriers are paramagnetic. 

It has been shown^ that two localized spins in a semiconductor can interact via the 
virtual carriers created in the (otherwise empty) bands by a laser of frequency smaller than 
the band gap. The effective interaction is a ferromagnetic Heisenberg coupling which we 
shall call optical RKKY (ORKKY). It is our contention that the ORKKY interaction can 
drive a diluted magnetic semiconductor into a ferromagnetic phase at low temperatures. In 
this paper, we develop a microscopic theory for this new type of ferromagnetism and we 
explore its properties. 

When the material is illuminated by a laser with electric field Ea, and frequency ul below 
the semiconductor band gap. Eg, a reactive dipolar energy is stored in the semiconductor: 

S = -J2xU^L)E^Ep (2) 

where x'aisi^L) is the real part of the retarded optical response function. In this case, as 
in many others where direct absorption processes are negligible, the laser field acts like a 
new tunable thermodynamic variable which changes the properties of the system. Xa/si^) 
depends on the inter-band transition energies which, due to the exchange interaction (P), 
depend in turn on the collective magnetization A4. In an illuminated sample, S does depend 
on A^. As we show below, S is minimized when the Mn spins are fully polarized and the 
system is ferromagnetic. Importantly, because Jtwl < Eg, the dissipative part x'Lisi^L) of 
the response function is zero, so that no real electron-hole pairs are created, and heating is 
strongly suppressed. 

We now derive S starting from the microscopic Hamiltonian. The relevant degrees of 
freedom are the electrons in the conduction band, holes in the valence band and the Mn spins. 
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The exchange interaction (P) is included in the VCMFA, which results in a spin splitting 
of the bands when the Mn spins are polarized^. For simplicity, we have neglected the 
important spin-orbit interaction influence on the valence bands; although these interactions 
do not change the qualitative physics, they do introduce a characteristic dependence on the 
laser light polarization that we do not discuss here. A classical laser field couples the valence 
and conduction bands. We only consider here the case of unpolarized laser light, where the 
two optically active (OA) interband transitions are driven with equal strength. We use the 
selection rule for the heavy holes at = so that an electron with spin a is promoted 
from the valence band to the conduction band. The resultant magnetization comes from 
the cooperative effects of the Mn spins through their interactions with the carriers, unlike 
optical orientation effects that can be induced by polarized light in these materials. 
The Hamiltonian reads H = Hc + + Hl + Vc with: 

ct Cr 



cr,k 



= 7lr;J2^{mm-Q) (3) 



2V 

where we have chosen the spin quantization axis along the collective magnetization axis 
Ai. The conduction and valence band dispersions are e| and respectively. The strength 
of the light matter coupling is quantified by the Rabi energy, Q = dcvEo where Eq is the 
amplitude of the electric field of the laser. The last term accounts for the Coulomb interaction 
in the dominant term approximation with p{q) = Peio) — Ph{(l), Pe{<l) = Xlfc o- '^jr^'^fc+q o- 
Phio) = u'^\ ^'^k+qa^ ^ i^l) = 7^ ^"^^ ^ stauds for the volume of the sample. 

The dynamics of the carrier density matrix under the influence of Hamiltonian Q in the 
spinless case with Je^h = has been studied previously by a number of authors, who treat 
Coulomb interactions in the Hartree Fock (HF) approximation^iiSiii. Numerical solution of 
the dynamical equationsii show that the density matrix reaches a steady state due to the 
inhomogeneous broadening provided by the different k states. Spontaneous emission and 
other phenomena will also contribute to the damping of the density matrix dynamics in a 
time scale T2 which is in the picosecond range. The HF steady state can be described as the 



vacuum of a dressed Hamiltonian-W and resembles very much a BCS state for excitons, where 
Q and the electron hole (eh) interaction act as the pairing force and 6q = Eg — ujl plays the 
role of the chemical potential. When 60 is positive and sufficiently large, the average density 
of band carriers vanishes identically when the laser is switched off. In that limit, the wave 
function of the ground state of the dressed Hamiltonian can be mapped to a coherent state of 
non- interacting bosonic excitons^ with binding energy e^, Bohr radius and ground state 
wave function ipis{r) corresponding to the electron hole Coulomb potential V{q) above. A 
crossover to the incoherent region, in which the density of carriers can be non zero when the 
laser is switched off occurs when 60 decreases towards smaller or even negative values. 

When the spin degree of freedom is included, the equations for the two optically active 
channels decouple in the HF approximation^ for an unpolarized laser. The two sets of 
equations for the carrier density matrix of the two OA channels, + and — , are identical 
to those of the spinless case2»i2iii except for the detuning, which acquires a spin dependent 
magnetic shift. It proves useful to define a detuning renormalized by the excitonic effects, 
6 = 6q — ex and a spin dependent detuning 



with J'(A^) = {Jh — Je)cMn2'- As a result of the magnetization, the detuning in one channel 
increases whereas in the other it decreases. For Ai = the coherent region (no absorption) 
occurs for 6 > C^^. In order to keep the semiconductor transparent in the ferromagnetic 
phase, we need 6 > J^{A4). For the DMS considered below, the magnetic shift is large 
enough so that the system is in the non-interacting exciton situation for which analytical 
results are available^.. For instance, the density of virtual excitons in each channel reads: 



This result is valid provided that n„a^^ < 0.1, which is the case for all the results presented 
below. This condition imposes both a lower bound for 6, which must be larger than ^min = 
\J^{A4 = S)\ = Jmax, and an upper bound for Q, which must be smaller than Quiax = 
1.12 X ((5 — 5inin)- The former is the condition for the material to remain non-absorbing in 
the ferromagnetic phase, whereas the latter is the limit for the validity of the independent 
exciton approximation. In this limit the reactive dipolar energy stored in the semiconductor 



(M) = Eg-UL-e,±J{M) = 6±J (M) 



(4) 



1 

4^ (5±)2 



(5) 
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because of the interaction with the laser reads: 



1 1 



(6) 
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which should be compared with Formally eq. ® is the ground state energy of the 
dressed Hamiltonian, in the non-interacting exciton limit. To explore the collective prop- 
erties of the Mn magnetic moments under the influence of the laser, we assume that the 
magnetic degrees of freedom are in thermal equilibrium and that their internal energy S is 
given by Eq. (jH]). The magnetization is determined by the competition between the dipo- 
lar energy density = ^ ^(-^) the entropy density in Qs = —CMnkbTS{A4), where 
S{M.) is the entropy per spin as a function of the magnetization^. When the paramagnetic 
to ferromagnetic phase transition is continuous, we can expand G [Ai) around = to 
obtain an analytical formula for the Curie Temperature: 

,,ro ^ '-^U. - AfcJ^^^. (7) 

We note that the same equation can be derived from the mean field theory using the ef- 
fective Heisenberg interaction defined by the optical exchange coupling derived in ref. 0. 
From equation ((7j) and the numerical evaluation of Tq shown below, we conclude that the 
important factors to achieve sizable Curie Temperature are compact excitons (small a^) 
and large Rabi energy Q. From this point of view, (Zn,Mn)S and (Zn,Mn)Se seem to be the 
most promising in the II- VI family, due to their small Bohr radii. 

In figure 1 we plot Qsi^Ai) and Gsi-^) for bulk Zno.988Mno.012S {xcs = 0.01) under the 
influence of a laser with Rabi energy = 5 meV and a detuning 6 = 41 meV— . For this 
concentration, the magnetic shift for the fully polarized case is Jmax =21 meV. The upper 
panel of figure 1-a shows Gsi-M.), which favours magnetic disorder {A4=0). In the middle 
panel of figure 1-a we see how the dipolar energy is minimized for a maximal magnetic order 
M = 5/2. In figure 1-b we plot the total Landau functional G (Ai) = + Qg, for a set 
of temperatures between 105 and 115 mK. For each temperature, the order parameter Ai 
is obtained by minimization of G {Ai) (magenta points in figure lb). The curve A4{T) so 
generated is shown in the lower panel of figure 1-a. A discontinuous phase transition from 
paramagnetic {At = 0) to ferromagnetic {At 7^ 0) occurs at Tc = 114 mK. In that panel we 
also show Ai{T) for some other detunings, 6 = 26 meV {Tc =778 mK), 6 = 71 meV {Tc =22 
mK). We see that Tc is a decreasing function of 6 and that the shape of the Ai{T/Tc) curve 



varies: our model predicts continuous and discontinuous phase transitions in different regions 
of the phase diagram. Conveniently, equation ((Tj) remains a good approximation for Tc even 
when the transition is discontinuous. 

The Curie temperature as a function of 6 and Q is shown in figure 2-A for Zno.988Mno.012S. 
From figure 2B it is apparent that ferromagnetism can be switched on and off, at a given 
temperature, by adjusting the intensity or the frequency of the laser. For a laser field 
with Rabi energy = 10 meV and 6 = 30.5 meV our theory predicts Tc = 1.5 Kelvin. 
Continuous wave laser excitation with Rabi splitting Q = 0.2 meV has been achieved^^ 
which would yield Tc ~ 50 mK for bulk Zno.988Mno.012S and 5 = 21.6 meV^. In the opposite 
limit, values of Q as large as 1.5 eV can be achieved with 5 femtosecond laser pulses^^. 
However, in order to observe ferromagnetic order the laser pulse must be longer than the 
transition switching time, which should presumably have the same order of magnitude as Ti, 
the Mn longitudinal relaxation time. In the absence of laser driven interactions, this quantity 
depends strongly on x, the Mn concentration, indicating that the superexchange interactions 
are the dominant relaxation mechanisnJ^. For x ^ 0.01, the laser driven interaction is 
dominant and presumably as efficient as the superexchange interactions in samples with 
larger Mn concentrations, where Ti can be of the order of 0.2 ns^^lS,. In contrast, the 
ferromagnetic interactions will switch on and off as fast as the laser pulse, on a femtosecond 
time scale with state of the art laser techniques. 

The full phase diagram in T and 6 would have six distinct regimes (see figure 2C). For 
S > Jma.x the semiconductor is always transparent and the system can undergo the purely 
coherent ferromagnetic phase transition at a Tc close to the estimates plotted in figure 
2A. For 0< 6 < Jmax, the system is transparent and paramagnetic at high temperature, 
but has the potential to become ferromagnetic and absorbing at low temperatures. In 
this situation, the system might return to the paramagnetic phase, because of heating, 
or might be ferromagnetic (due to the standard RKKY interaction mediated by the real 
photocarriers). Finally, for negative detuning the system is always absorbing, and photo- 
induced carriers can lead to ferromagnetism through the same mechanisms that apply for 
equilibrium carriers, leading to a transition at ksTc = ^^^^^^ f^^^^l^iij-i) {Jsd^c + Jpcff^v) 
where n is the density of photocarriers. Ferromagnetism induced by (real) photocarriers has 
been observed in (Cd,Mn)Te heterostructures at 2 Kelvini^ 

In summary, we propose that a diluted magnetic semiconductor exposed to a laser below 
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the absorption threshold will undergo a ferromagnetic phase transition at low temperatures. 
The ferromagnetism takes place because the reactive dipolar energy induced in the semicon- 
ductor is minimized by the ferromagnetic configuration. Microscopically, the spins interact 
with each other via an optical RKKY interaction mediated by virtual carriers®. The ferro- 
magnetic interactions are switched on and off as fast as the laser. The mechanism of this 
laser induced ferromagnetism relies on interband optical coherence and is radically different 
from other types of carrier mediated ferromagnetism. 
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CAPTION Figure 1 

Energy contributions as a function of A4. Figure 1-A. Upper panel: Qsi-M.) for two 
different temperatures. Middle panel: ^^y^ for 5 = 41 meV and = 5 meV. Lower panel: 
A4(T/Tc) for Q = 5 meV and three different values of 6. Figure 1-B: total Landau energy 
Q{Ai) for the Mn, for a set of temperatures nearby Tc for 6 = 41 meV and f2 = 5 meV. 

CAPTION Figure 2 

(A). Contour map for transition temperature as a function of the Rabi energy Q and the 
detuning 6, for (Zno.99,Mno.oi)S. White region: outside validity range of linear response. (B) 
A4{fl,6) at a fixed temperature, kBT=0.5 Kelvin for Zno.988Mno.012S. (C) Schematic {ksT, 
S/Jma.x) phase diagram, for fixed laser intensity. In the left region {6 < 0) (right region 
{6 > Jmax) ) the system is always absorbing (coherent). The red line in the coherent side 
signals the ferromagnetic instability described in this work. The red line in the absorbing 
part corresponds to photocarrier mediated ferromagnetism T,r*i^. The middle region shows a 
possible transition from a paramagnetic-coherent towards a ferromagnetic-absorbing phase. 



10 




Figure 1 



11 



A 




30 35 40 45 50 55 60 

8 (meV) 




